Antiferromagnetic Heisenberg ladders in a staggered magnetic field 
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We study the low-energy excitations of the spin- 1/2 antiferromagnetic Heisenberg chain and A'^-leg 
{N=2, 3, 4) ladders in a staggered magnetic field hs- We show that kg induces gap and midgap states 
in all the cases and examine their field scaling behavior. A modified boundary scheme is devised to 
extract accurate bulk excitation behavior. The gap values converge rapidly as A'^ increases, leading 
to a field scaling exponent 7 = 1/2 for both the longitudinal and transverse gaps of the square 
lattice {N — > oo). The midgap states induced by the boundary edge effects share the bulk gap 
scaling exponents but their overall scaling behavior in the large-A limit needs further investigation. 
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Heisenberg spin ladders have attracted considerable in- 
terest for their fascinating properties due to strong quan- 
tum fluctuations and their unique structural character as 
a crossover platform from one to two dimensions^ It is 
known that spin ladders with odd and even numbered 
legs have drastically different behavior; the former have 
gapless spin excitation spectra while the latter have fi- 
nite spin gapsiSi^ This disparity complicates consider- 
ably attempts to extrapolate the results of ladders to two- 
dimensional systems. An important development in the 
study of low-dimensional spin systems is the observation 
of an unexpected magnetic field induced gap in the low- 
energy excitation spectrum of copper benzoate, a quasi- 
one-dimensional spin chain material.'^ Recently, the field 
induced gap has been reported in more spin chain ma- 
terials such as Yb4As3,^'^'* CuCl22(dimethylsulfoxide)^, 
and [Pyrimidine-Cu(N03)2(H20 )2]„(CuPM).i'' 

A gap in the low-energy excitation spectrum of a spin- 
1/2 Heisenberg chain can be induced by an effective stag- 
gered magnetic field >ii The staggered field may orig- 
inate from the staggered gyromagnetic tensor or/and 
the Dzyaloshinskii-Moriya (DM) interaction when an ex- 
ternal magnetic field is appliedii^ Using a local gauge 
transformation and neglecting small anisotropic terms in 
both the Heisenberg exchange and the Zeemann splitting 
terms, an effective Hamiltonian for spin chain materials 
with the DM interaction can be written asl?> 
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J^i ■ Si+i — HS^ — hs{—iyS![ 
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where H and hs are the uniform and staggered magnetic 
field, respectively. This effective Hamiltonian has been 
mapped onto the sine-Gordon model using the bosoniza- 
tion techniquc^^ to obtain an analytic form for the spin 
gap as a function of the magnetic field. Furthermore, 
it was shown^^ that the uniform field does not change 
the qualitative scaling behavior induced by the staggered 
field and the scaling function derived for H=0 provides a 
good explanation for the experimental results at the low- 



field region and its validity is supported by the density 
matrix renormalization group (DMRG) calculationsiiSii^ 
However, the effective field theory is only applicable when 
the field is not too strong. The recent DMRG calcula- 
tions unveiled different field-dependence at higher fields 
and a crossover in the intermediate field regime^^*i^ 

Additional intriguing phenomena induced by the stag- 
gered field have been recently reported for quasi-one- 
dimensional spin chain materials BaCu2Si20'r^ and 
CuCl2 ■ 2{dimethylsulf oxide) (CDC).^ The neutron 
scattering measurement on CDC and electron spin reso- 
nance measurement on BaCu2Si207 show that the field 
dependence of the induced gap deviates from the sine- 
Gordon model prediction for the spin chains. It suggests 
that the interchain interaction may play an important 
role. Studies on the effect of the interchain interaction 
is therefore needed to clarify the fundamental physics 
and to establish the nature of the low-energy excitation 
in these new materials. Moreover, there is also consider- 
able interest in searching for a reliable extrapolation from 
multi-leg ladders to the square lattice-iSJI This should 
be achievable for spin chain materials with the DM in- 
teraction since, unlike the hs=0 case, ladders with both 
odd and even numbered legs show field induced gaps with 
similar scaling behavior under the staggered field. 

In the present work, we study spin-1/2 antiferromag- 
netic Heisenberg ladders with the following Hamiltonian: 
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H — </Sa,i • Sa,j+1 + Jl_^a,i ' ^a+l,i 

a—1 i=l a—1 i—1 

N L 

+ EE(-l)^^"^«^a,.. (2) 



where N (L) is the number of legs (rungs). We consider 
the case of isotropic coupling, i.e., J — J± = 1 and em- 
ploy the DMRG method^^'^^''^'^ to study its low-energy 
properties. We kept up to 800 states for ladders with 
up to 300 rungs in our computations. The truncation 
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errors are less than 10~® in all the cases. Although the 
staggered field breaks the SO (3) symmetry, the total 
remains conserved. The spin gaps for the longitudinal 
and transverse branches, A£ and A7-, are defined as: 

l^c{L) = E^{L,Q) - E^{L,Q) (3) 
Ar(i) =-Bo(i,l)-£^o(i,0) (4) 

where Eq{L, S^) and Ei(L, S^) are the ground-state and 
first excitation energy in the sector. In DMRG calcu- 
lations, numerical accuracy is usually much higher with 
the use of the open boundary condition (OBC) than 
that with the periodic boundary condition (PBC). Re- 
cent studies^ ^ reveal that there generally exist midgap 
states induced by the edge excitations in open-end spin 
chains in the presence of a staggered magnetic field. Our 
calculations show that midgap states also exist in spin 
ladders (see below). This complicates the process of ex- 
tracting the bulk excitations. To address this issue, we 
devised a modified boundary scheme (MBS) to first iden- 
tify and then move the midgap states away from the low- 
energy spectrum. It ensures a reliable extraction of the 
bulk excitation gap while allowing an accurate descrip- 
tion of the midgap states. 

The idea behind the MBS is to push the edge excita- 
tions up in energy and leave only the bulk excitations in 
the low-energy spectrum. This is achieved by introduc- 
ing an edge parameter in the Hamiltonian that systemat- 
ically drives up the edge excitation. Similar treatments 
have been proposed in previous DMRG calculations with 
OBCi^2i2i24 Here we introduce a continuous monotonic 
function /(a;),0 < a; < 1, where /(O) > 1 and /(I) = 1, 
to rescale the Hamiltonian parameters 

{a, l<i < M 

1, M <i<L- M 

CL+i-i, L-M + l<i<L 

where Y represents J, hg, or J± and c™ — f{m/M), 1 < 
m < M where M is the number of the edge sites on which 
the parameters are adjusted. We choose the functional 
form f{x) = 1 -|- a(l -|- cos(7rx)) where a > is the only 
adjustable parameter (a = corresponds to OBC). This 
choice is certainly not unique, but it satisfies the require- 
ment for simplicity and effectiveness in removing edge 
excitations from the low-energy spectrum. We have run 
extensive tests and found this scheme work very well for 
ladders with both ferromagnetic and antiferromagnetic 
couplings at all field strengths. 

We demonstrate the implementation of MBS in the 
case of the two-leg ladder. Figure ^ shows the gap and 
midgap states versus 1/L with various a at /is = 0.6. 
The gap obtained under PBC represents the bulk value. 
A midgap state appears in the open-end (OBC) case 
due to topological edge effects i^^iSSiSl As a increases, the 
midgap state is gradually pushed up in energy and even- 
tually enters the continuum spectrum, leaving the bulk 
gap state as the lowest excitation in the large-L limit. 
We performed the same calculations on the three- and 
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FIG. 1: The longitudinal (top) and transverse (bottom) gap 
and midgap states versus 1/L at various ct values for the two- 
leg ladder at hs=Q.& with M = 10. 

four-leg ladders and found that midgap states exist for 
all the cases. With MBS one can clearly identify these 
bound states inside the field induced gap for a detailed 
study and then systematically push them up in energy 
and extract the bulk excitations with proper values of a. 
It is noted that a large a is needed for cases where the 
edge state lies well below the bulk gap. In some special 
cases where the edge states are degenerate in the large-L 
limit with the ground state as in spin-1 and 2 Heisenberg 
chains, special measurea^^*^ need to be taken to elimi- 
nate the degeneracy before applying MBS. In the calcu- 
lations presented below we use a = 2.0 with M — 10 to 
ensure good convergence in all the cases. 
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FIG. 2: The longitudinal and transverse gap versus the stag- 
gered magnetic field hs for the spin-1/2 chain (A^=l) and 
A^-leg (A^ = 2, 3, 4) ladders. The fitted scaling curves are also 
shown for = 1,3,4 as the dashed, dotted and solid lines, 
respectively. The insets illustrate the gaps for hs £ [0,0.1]. 



We now turn to the study of the staggered field induced 
gaps in the spin chain and ladders and the extrapolation 
toward the square lattice. Figure [3 shows the longitudi- 
nal and transverse gaps for the chain (7V=1) and iV-leg 
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(7V=2, 3, 4) ladders. When hs = the spin excitations 
for N—1 and 3 cases are gapless. For the two- and four- 
leg ladders, the zero (staggered) field gaps are 0.502J 
and 0.165 J, respectively. These are in agreement with 
previous DMRG results^ and Monte Carlo dataJl For 
the four-leg ladder, the present gap value is expected to 
be much more accurate than those previously reported 
since we have kept a very large number (up to 800) of 
states in the DMRG calculations. In the weak field limit, 
the gaps for the ladders with even-numbered legs show a 
different field scaling behavior from that for the ladders 
with odd- numbered legs (see the insets in FigO. In par- 
ticular, it has been showni2iiii4 that the field induced 
gap for the spin- 1/2 chain scales as hj with a scaling ex- 
ponent 7 = 2/3 and an additional logarithmic correction 
for small h^. Meanwhile, the gaps converge rapidly be- 
yond the weak field limit as N increases as seen in Fig. 

m 

In the large- iV limit, the gaps scale to zero at hg — 
as expected for the two dimensional case. It is also antic- 
ipated that the gaps of the ladders can be extrapolated 
to the two-dimensional case beyond the weak-field region 
given their smooth change and rapid convergence with 
increasing TV. To this end, we fit the field dependence of 
the gaps using the analytical scaling function 

A^Ao + fi{h,-h,,y, (5) 

for both transversai^ and longitudinal branchiSSiSfliSi As 
expected, Aq and hgg are zero for odd N within the nu- 
merical accuracy of the computation and scale toward 
zero monotonically and rapidly for even N. Fitting to 
the DMRG results indicates that hgg is already vanish- 
ingly small for N=A. In Table I, we show the direct fitting 
data for N = 1,3,4 using Eq. and the fitting curves 
are shown in Fig. [3 Meanwhile, due to the large gap at 
hs ~ 0, a, good fitting using Eq. (0) cannot be achieved 
for N—2; instead, the scaling parameters for N=2 are 
obtained from the extrapolation curves obtained from fit- 
ting the results for N=l, 3, 4 (see below). 



TABLE I: Fitting coefficients and exponents for A'^ — 1, 2, 3, 4 
and oo. Those for the two-leg ladder are taken from Fig. |3 
at A'^ = 2 and P and 7 for A = 00 from Fig. |3(see text). 
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1 


2 


3 


4 


00 


p¥ 

7? 


2.97 
0.68 
1.88 
0.61 


3.84 
0.56 
2.10 
0.56 


3.98 
0.53 
2.16 
0.54 


4.02 
0.52 
2.19 
0.53 


4.03 ±0.09 
0.50 ±0.01 
2.27 ±0.01 
0.50 ±0.01 



We extrapolate the longitudinal and transverse gaps 
7^^ r ^^'^ Pc T toward the large- A^ (two dimensional) 
limit using the second order polynomial fitting: a{N) — 
ao + aiA^^-'^ + a2N^'^, where a denotes either 
13^ T given in Table I with A^ = 1,3,4. This fitting is 
illustrated in Fig. Oand the extrapolated ^^q- and 



are given in Table I. From this procedure, we obtain the 
longitudinal and transverse gaps for the isotropic anti- 
ferromagnetic Heisenberg square lattice with a staggered 
magnetic field 

A£ = (4.03 ± 0.09)/i°-^°±°-°\ (6) 
Ar = (2.27±0.01)/i°-^°=^°°^ 

Our numerical results suggest a common exponent 7 for 
both branches. It is also interesting to note that a recent 
field-theoretical study^^ of coupled spin-1/2 antiferro- 
magnetic Heisenberg chains with a leg-independent stag- 
gered field also shows that the field induced gap scales 
with the staggered magnetic field as hs . Possible con- 
nections and general implications of the common scaling 
behavior in these different spin lattice models deserve 
further analytical investigation. 
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FIG. 3: The extrapolation of the scaling exponents and coef- 
ficients of the field induced gaps from A-leg ladders toward 
the square lattice {N 00). 



Finally, we examine the scaling behavior of the midgap 
states in spin ladders with a staggered magnetic field. 
As shown in Fig. 4, midgap states exist in all the cases 
studied and the scaling exponents 7^ and 7^ for the 
bulk gaps provide a good fit for the midgap states. How- 
ever, the scaling coefficients fi^^ and (3^ show an oscil- 
latory pattern with increasing A^ (up to A^ = 4 studied 
here). It indicates that these coefficients have not prop- 
erly converged at A^ = 4, probably due to a different 
(compared to the bulk excitations) energy scale associ- 
ated with the boundary edge excitationsiSSiSSiSl Conse- 
quently, the question of whether the midgap states would 
persist in the large- A^ limit remains open at present. 
Meanwhile, we observe that although the bulk low-energy 
behaviors of the spin-1/2 two-leg ladder are similar to 
those of a spin-1 chain,"^^ their boundary edge excita- 
tions and the field scaling behaviors are quite different!^ 
The spin-1 chain has a four-fold degenerate ground state 
that splits in a staggered field and turns (partly) into the 
lowest midgap states.-^ The two-leg ladder has a non- 
degenerate ground state and its transverse midgap states 
come from the boundary edge excitations that originally 
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FIG. 4: The scaling fit for tlie transverse and longitudinal 
gaps and the midgap states for the ladders (top to bottom 
panels: A^=2— 4). The filled and open symbols represent the 
gap and midgap states, respectively. The exponents and co- 
efficients for 1^=2 are taken from Fig. |3 



overlap with the bulk excitation continuum. However, 
when one adds other interactions to the standard two- 
leg ladderiSiS, the resulting low-energy properties, ie. gap 



and midgaps, in the Haldane phase can be the same as 
the S" = 1 chaii^. 

In summary, we have carried out a systematic study 
of the low-energy excitations of spin-1/2 antiferromag- 
netic Heisenberg chain and ladders in a staggered mag- 
netic field and obtained the field scaling behavior of both 
longitudinal and transverse gaps. A modified boundary 
scheme has been devised to extract the bulk excitation 
behavior using the more accurate open boundary option 
in the DMRG calculations. The calculated bulk gaps con- 
verge rapidly with increasing number of legs above weak 
field regime; it allows a reliable extrapolation to obtain 
the field scaling behavior of the gaps for the isotropic 
square lattice. We also examined the midgap states in 
the spin ladders induced by the topological edge effect. 
The midgap states share the same scaling exponents with 
the bulk gaps but the scaling coefficients show an oscilla- 
tory pattern with increasing number of legs (up to iV=4) . 
As a result, the overall scaling behavior of the midgap 
states in the large- A'^ limit needs further investigation. 
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